ABSTRACT
BACKGROUND: To identify host genetic variants (SNPs) associated with COVID-19 disease severity, a number of genome-wide association studies (GWAS) have been conducted. Since most of the identified variants are located at non-coding regions, such variants are presumed to affect the expression of neighbouring genes, thereby influencing COVID-19 disease severity. However, it remains largely unknown which genes are influenced by such COVID-19 GWAS loci. METHODS: CRISPRi (interference)-mediated gene expression analysis was performed to identify genes functionally regulated by COVID-19 GWAS loci by targeting regions near the loci (SNPs) in lung epithelial cell lines. The expression of CRISPRi-identified genes was investigated using COVID-19-contracted human and monkey lung single-nucleus/cell (sn/sc) RNA-seq datasets. FINDINGS: CRISPRi analysis indicated that a region near rs11385942 at chromosome 3p21.31 (locus of highest significance with COVID-19 disease severity at intron 5 of LZTFL1) significantly affected the expression of LZTFL1 (P<0.05), an airway cilia regulator. A region near rs74956615 at chromosome 19p13.2 (locus located at the 3' untranslated exonic region of RAVER1), which is associated with critical illness in COVID-19, affected the expression of RAVER1 (P<0.05), a coactivator of MDA5 (IFIH1), which induces antiviral response genes, including ICAM1. The sn/scRNA-seq datasets indicated that the MDA5/RAVER1-ICAM1 pathway was activated in lung epithelial cells of COVID-19-resistant monkeys but not those of COVID-19-succumbed humans. INTERPRETATION: Patients with risk alleles of rs11385942 and rs74956615 may be susceptible to critical illness in COVID-19 in part through weakened airway viral clearance via LZTFL1-mediated ciliogenesis and diminished antiviral immune response via the MDA5/RAVER1 pathway, respectively. FUNDING: NIH.
Subject(s)
COVID-19/genetics , CRISPR-Cas Systems , Genetic Loci , Polymorphism, Single Nucleotide , Ribonucleoproteins/genetics , SARS-CoV-2/genetics , Transcription Factors/genetics , Animals , COVID-19/metabolism , Chromosomes, Human, Pair 19/genetics , Chromosomes, Human, Pair 19/metabolism , Chromosomes, Human, Pair 3/genetics , Chromosomes, Human, Pair 3/metabolism , Databases, Nucleic Acid , Genome-Wide Association Study , Haplorhini , Humans , RNA-Seq , Ribonucleoproteins/metabolism , SARS-CoV-2/metabolism , Transcription Factors/metabolismABSTRACT
Genome-wide association studies have identified lung disease-associated loci; however, the functions of such loci are not well understood in part because the majority of such loci are located at non-coding regions. Hi-C, ChIP-seq and eQTL data predict potential roles (e.g. enhancer) of such loci; however, they do not elucidate the molecular function. To determine whether these loci function as gene-regulatory regions, CRISPR interference (CRISPRi; CRISPR/dCas9-KRAB) has been recently used. Here, we applied CRISPRi along with Hi-C, ChIP-seq and eQTL to determine the functional roles of loci established as highly associated with asthma, cystic fibrosis (CF), chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF). Notably, Hi-C, ChIP-seq and eQTL predicted that non-coding regions located at chromosome 19q13 or chromosome 17q21 harboring single-nucleotide polymorphisms (SNPs) linked to asthma/CF/COPD and chromosome 11p15 harboring an SNP linked to IPF interact with nearby genes and function as enhancers; however, CRISPRi indicated that the regions with rs1800469, rs2241712, rs12603332 and rs35705950, but not others, regulate the expression of nearby genes (single or multiple genes). These data indicate that CRISPRi is useful to precisely determine the roles of non-coding regions harboring lung disease-associated loci as to whether they function as gene-regulatory regions at a genomic level.